case = "mz.1440x720x30"

f = addfile(case + ".h0.nc", "r")

itime = 25

wks = gsn_open_wks("pdf", case)

res = True
res@gsnDraw = False
res@gsnFrame = False
res@cnLinesOn = False
res@cnFillOn = True
res@cnFillPalette = "amwg256"
res@mpOutlineOn = False
res@mpCenterLonF = 180

u700 = linint1_n_Wrap(f->ph(itime,:,:,:), f->u(itime,:,:,:), False, (/70000.0/), 0, 0)
u700&LEV = (/u700&LEV / 100/)
u700@long_name = u700@long_name + "@700hPa"

v700 = linint1_n_Wrap(f->ph(itime,:,:,:), f->v(itime,:,:,:), False, (/70000.0/), 0, 0)
v700&LEV = (/v700&LEV / 100/)
v700@long_name = v700@long_name + "@700hPa"

t700 = linint1_n_Wrap(f->ph(itime,:,:,:), f->t(itime,:,:,:), False, (/70000.0/), 0, 0)
t700&LEV = (/t700&LEV / 100/)
t700@long_name = t700@long_name + "@700hPa"

plot = new(3, graphic)
plot(0) = gsn_csm_contour_map(wks, u700(0,:,:), res)
plot(1) = gsn_csm_contour_map(wks, v700(0,:,:), res)
plot(2) = gsn_csm_contour_map(wks, t700(0,:,:), res)

res := True
gsn_panel(wks, plot, (/1,3/), res)
